Dilepton production in a schematic causal viscous hydrodynamics 
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Assuming that in the hot dense matter produced in relativistic heavy-ion collisions, the energy 
density, entropy density, and pressure as well as the azimuthal and space-time rapidity components 
of the shear tensor are uniform in the direction transversal to the reaction plane, we derive a set 
of schematic equations from the Isreal-Stewart causal viscous hydrodynamics. These equations are 
then used to describe the evolution dynamics of relativistic heavy-ion collisions by taking the shear 
viscosity to entropy density ratio of l/47r for the initial quark-gluon plasma (QGP) phase and of 
ten times this value for the later hadron-gas (HG) phase. Using the production rate evaluated 
with particle distributions that take into account the viscous effect, we study dilepton production 
in central heavy-ion collisions. Compared with results from the ideal hydrodynamics, we find that 
although the dilepton invariant mass spectra from the two approaches are similar, the transverse 
momentum spectra are significantly enhanced at high transverse momenta by the viscous effect. 
We also study the transverse momentum dependence of dileptons produced from QGP for a fixed 
transverse mass which is essentially absent in the ideal hydrodynamics, and find that this so-called 
transverse mass scaling is violated in the viscous hydrodynamics, particularly at high transverse 
momenta. 
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I. INTRODUCTION 

The ideal hydrodynamics without shear viscosity has 
been quite successful in describing the anisotropic flow of 
particles in heavy-ion collisions at the Relativistic Heavy 
Ion Collider (RHIC) Since the viscosity is inversely 

proportional to the scattering cross section between con- 
stituent particles, the applicability of the ideal hydro- 
dynamics at RHIC provides a strong evidence that the 
quark-gluon plasma formed at RHIC is a strongly coupled 
one (sQGP). The study based on the Ads/CFT gauge- 
gravity duality has, however, suggested that the shear 
viscosity to entropy density in the QGP cannot be smaller 
than l/47r. Small viscosities of QGP have also been ob- 
tained in studies based on either the quasi-particle model 
that fits the equation of state from the lattice gauge cal- 
culations Q or the pQCD including both gluon clastic 
and radiative scatterings that gives a good description of 
measured elliptic flows at RHIC On the other hand, 
the viscosity of hadronic matter has been found to be 
much larger in theoretical studies [^-Qj about an order 
of magnitude larger than the lower bound predicted by 
the Ads/CFT gauge-gravity duality. 

Including a small viscosity in the hydrodynamics has 
led to an improved description of measured anisotropic 
flows of hadrons at large transverse momenta A 
nonzero viscosity also affects particle momentum distri- 
butions in the hot dense matter produced in heavy-ion 
collisions, resulting in a deviation from thermal equilib- 
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rium during its expansion, and this is particularly so for 
particles of high transverse momenta. Moreover, the vis- 
cosity can change the evolution dynamics of produced 
hot matter in relativistic heavy-ion collisions. Because of 
the heat generated by viscosity, cooling of the hot mat- 
ter becomes slower, leading to a slower decrease of its 
temperature and thus a larger transversal but a slower 
longitudinal expansion compared to the case of the ideal 
hydrodynamics [l0| . 

The viscous effect can further affect dilepton produc- 
tion in relativistic heavy-ion collisions, which has been 
suggested as a possible tool to probe the properties of 
quark-gluon plasma 11 - [T6j as well as those of hot dense 
hadronic matter [i^^ill. As shown in Ref. [23 |. the vis- 
cosity modifies significantly the transverse momentum 
spectrum of dileptons produced in relativistic heavy-ion 
collisions, although not much the invariance mass spec- 
trum. This study was based on the non-causal Navier- 
Stokes viscous hydrodynamics, which is known to have 
instabilities in numerical simulations [25j . and only in- 
cluded the viscous effect on dilepton production due to 
the leading order correction from modified particle distri- 
butions. In the present paper, we extend the study by us- 
ing the causal Israel-Stewart viscous hydrodynamics and 
including also the effect of the second-order correction 
from modified particle distributions on dilepton produc- 
tion. To simplify calculations, we assume that in the hot 
dense matter produced in relativistic heavy-ion collisions, 
the energy density, entropy density, and pressure as well 
as the azimuthal and space-time rapidity components of 
the shear tensor are uniform in the direction transversal 
to the reaction plane and derive a set of schematic equa- 
tions as in Refs. [2^, [13] for the ideal hydrodynamics. 

This paper is organized as follows: In Sec. [Hi general 
causal viscous hydrodynamic equations are introduced. 
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These equations are simplified in Sec. IIIII by assuming 
that all thermal quantities as well as the azimuthal and 
space-time rapidity components of the shear tensor are 
uniform in the transversal direction. In Sec. IIVI the 
quasi-particle model of massive quarks and gluons for the 
quark-gluon plasma and the resonance gas model for the 
hadronic matter are introduced for describing the equa- 
tion of state of the hot dense matter produced in relativis- 
tic heavy-ion collisions. The schematic viscous hydrody- 
namic equations are then solved in Sec. |V]for heavy ion 
collisions at RHIC energies to show the viscous effect on 
the evolution dynamics and the distribution of particle 
transverse momentum spectra. In Sec. IVIl dilepton pro- 
duction is studied in the viscous hydrodynamics to find 
out the effect of viscosity on the dilepton invariant mass 
and transverse momentum spectra. Finally, discussions 
and summary are given in Sec. lVIll Details on the deriva- 
tion of the dilepton production rate using the modified 
particle distributions from the viscous hydrodynamics is 
given in the Appendix. 

II. THE VISCOUS HYDRODYNAMICS 

In hydrodynamic description of relativistic heavy-ion 
collisions, the hot dense matter is characterized by its 
net charge currents and energy-momentum tensor. Since 
particles at midrapidities are largely produced ones, their 
net charge currents are essentially zero and can be safely 
neglected. In the Landau and Lifshitz frame, which as- 
sumes that the four- vector velocity u'^ = 7(1, w) is par- 
allel to the energy flow and the heat conductivity is 
zero [28|. the energy- momentum tensor can be written 
as [231 



T""" = (e + p)u''u'' - pg"" + tt"" , 



(1) 



where e and p are the energy density and pressure, re- 
spectively, and 'K^'^ is the traceless symmetric shear ten- 
sor. At midrapidity particles follow essentially the boost- 
invariant expansion along the longitudinal direction [soj , 
i.e., the longitudinal flow velocity is equal to z/t, if it 
starts at z = t = 0. Furthermore, the transverse flow ve- 
locity is independent of the azimuthal angle 4> in central 
heavy-ion collisions. In the (r, r, (j), rf) coordinate system 
defined by 



7] 



1, t + Z 

— In , 

2 t-z 



r — \J + y^, = tan ^{y/x)^ 



(2) 



only T^^, T"', T^'', T'''' and T"^"^ components of the 
energy-momentum tensor, and tt'^'^, tt''^, tt'^'", tt'''' and 
T^'txt' of the shear tensor are non-zero in central heavy-ion 
collisions. For the energy-momentum tensor, they are 
given by 



T'-'- = {e + Pr)ul+Pr 



(3) 



where Pr = p — t'^tt^^ — r'^n'^'l' is the effective radial pres- 
sure. The azimuthal and space-time components of the 
shear tensor r'^n'^'t' and t'^tt^^ are the only independent 
ones as the others can be related to them according to 



VrTT 



2 rr 



(4) 



where the first two equations are derived from w^tt'"' = 
and the last one from the traceless property tt^ = 0. The 
shear tensor components tt"^"^ and tt'''' are boost-invariant 
in the radial direction and satisfy following simplified 
Israel-Stewart equations: [2^ 



{dr + Vrdr)7r"" = 
{dr + Vrdr)TT''"'' = - 

where 
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(5) 
(6) 



9 ^ d ■ U ~ —dr{TJr) H dr{rVrJr) 

T r 

with rjs and r^r being the shear viscosity and the relax- 
ation time for the particle distributions, respectively. 

From the energy-momentum conservation conditions 
dfiT^" = 0, we then obtain following viscous hydrody- 
namic equations for the produced fire-cylinder: 

-dArr-) + ldr.{rrn = --(p + t^tt'"'), (7) 
T r T 



1 



1 



1 



-{p + r^TT't"t'). 



(8) 



Furthermore, the condition u^{T'^J^) — 0, where the flow 
velocity ( ) = (7/ cosh 77, "fVr, 0, 0) reduces to 

(7r, 7rfr, 0, 0) with jr = 1 / ^^1 — in midrapidities, 
leads to 



T r 1 



(9) 



where s = (e +p)/r is the local entropy density in the 
hot dense matter. Eq. (jOj shows that a nonzero shear 
tensor affects the entropy density of the matter. 



III. A SCHEMATIC VISCOUS 
HYDRODYNAMICS 

If all thermal quantities like energy density, temper- 
ature, entropy density, and pressure as well as the az- 
imuthal and space-time components of the shear tensor 
are uniform along the transverse direction in the hot 
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dense matter produced in heavy-ion collisions, we can 
then simplify the causal viscous hydrodynamic equations 
by integrating them over the transverse area [261 [27l |. 
Specifically, we integrate Eqs. ([7]) and (O as well as 
Eqs. dU and ([6]) multiplied by and 7^?"^, respec- 

tively, over the transverse area. In terms of Trjj = t^tt'''' 
and ttJJ = t^tt'''', this leads to 



dr{AT{T^^)) = -{p + ^^^A, 



(10) 



T 



a,(Ars(7,)) 



A{lr) 



5r(A(7.)<') 



R 

dr{A{^r)) 

2?7, 



(11) 



W _ (tt) 

3 T 



(12) 



5.(^(7.) 7:l)-\dAA{lr)) + 2Al^-^ 



A 



7,.w,. 



(13) 



where A = 7ri?^ with R being the transverse radius of 
the uniform matter and (• ■ •) denotes average over the 
transverse area. Assuming that the radial flow velocity 
is a linear function of the radial distance from the center. 



i.e 

1 



r^^Vr = jRR{r/R), where R ~ dR/dr and 7_r 
1 — we then have 

(7.^)^1 + ^^^^ 
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IrVr 



37|i?2 



(7^ - 1) 



1R.R 

R ■ 



(14) 



Since the energy density e and pressure p are related by 
the equation of state of the matter through its temper- 
ature T, Eqs. (fT0|) -([T3 | are thus four simultaneous equa- 



tions for T, i?, TT? and tt'' 



IV. THE EQUATION OF STATE 

For the equation of state of QGP, we use the quasi- 
particle model of Ref. [3l|, which assumes that the QGP 
is composed of noninteractig massive quarks and gluons. 
In terms of the temperature T of QGP, their masses are 
given by 

Ni\g\T)T^ 
3 6 



m 



2 _ 9\T)T^ 



« 3 ' 

where the strong coupling constant g{T) is given by 

F{T,T,,A) = 



(15) 



(UA^, -27V/) In F2(T,T„ A)' 

18 T Tc 

18.4e(^/^c)V2 + i^X' 



with Tc = 170 MeV, TJA = 1.05, A^c = 3 and Nf = 
3. The pressure, energy density, and entropy density of 
QGP arc then given, respectively, by 

^'(^) = E/V / dkMT)--B{T) 
^ Gn^ Jo Ei 

^ Po{T) - B{T) 
<T) = E #T / dkPMT)E, + B{T) 



27r^r Jo 



, (16) 



with mi{T) and gi being, respectively, the thermal mass 
and degeneracy factor of parton species i. The parton 
distribution function is denoted by 



1 



(17) 



with the plus and minus signs in the denomi nator for 
quarks and gluons, respectively, and Ei = yj mf -\- k^. 
For the bag pressure B{T), it is determined from the 
relation s = dp/dT such that 



dpo dm'j{T) 



(T) dT ' 



(18) 



where Bq is the bag pressure at Tc and is taken to be 0.095 
times the energy density at this temperature in order to 
keep the pressure continuous at Tc. A similar value of 
Bq has been used in Ref. [32| for the case of nozcro but 
small baryon chemical potential. 

For the HG phase, we use the resonance gas model that 
includes both stable hadrons and their resonances up to 
1.5 GeV for mesons and 2.0 GeV for baryons. Its pres- 
sure, energy density and entropy density can be similarly 
evaluated as those in Eq. (|16p for the QGP, except that 
the bag constant is not present and the hadron masses 
are taken to be their values in free space. 

Because of the larger entropy density in QGP than 
in HG at Tc, a mixed phase of constant temperature Tc 
is introduced during the transition between these two 
phases of matter. In terms of the fraction / of HG in the 
mixed phase, the entropy of the mixed phase is 



(19) 



where and are, respectively, the entropy density of 
HG and QGP at Tc. Similar relations hold for the energy 
density and pressure. 
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We refer the reader to Ref. [33| for details of the above 
equations of state for QGP and HG. 

For the shear viscosity, we take its ratio with respect 
to the entropy density to be l/47r for QGP as given by 
the Ads/CFT gauge-gravity duahty Q and ten times this 
value for HG as determined from the hadronic transport 
model 01 . For the relaxation time in Eqs. ^ and (O, 
wc use the assumption ry/r^ = sT/3 for both QGP 
and HG. 



V. HEAVY-ION COLLISION DYNAMICS IN 
THE SCHEMATIC VISCOUS HYDRODYNAMICS 



To apply the schematic viscous hydrodynamic equa- 
tions, Eqs. ([TU1) - P^ . to heavy-ion collisions, wc divide 
the time into infinitesimal intervals. These equations 
then become 



(7^)„+ie„+i -t- (7^u^)„+i 



Pn+l ~ «,')"+! - (7r^)n+l 



ilr) n+1 



Tn+1 



At 



At 



TnSn+1 - (TTpn - (7r,'J)„ 



ilr) n+1 

RlAr 



(7r)r 



TnSn - { 1 



2(7fli?)„AT\ ^ 



{Tr)nRn j ^ 

2(7r)„ 1 



(7^i^)„Ar Ar- 

{lr)nR. - 



('^7r)n 



{lr)r 



{1t\ 



n+1 



fl|Ar 



IrVr 



(;^}<':).-lfe),.{(^-a<->"-(^ 



(20) 



(21) 



(22) 



(23) 



where Ar — t„+i — t,i with the subscript denoting the 
time step and (90/9r)„ = (On+i — On)/ At. 

Eqs. ([20)) - ((23|) are solved by using following relations 
between the (n + l)th and nth time steps in the energy 
density, pressure and entropy: 



where 
ds 
df 



T 



de 
'df 



dp 

Pn+l = Pn + ^ 

ds 

Sn+l — Sn + 



{Tn+1 — Tn), 

I 

{Tn+1 — Tn), 

I 

{Tn+1 ^ Tn), 



E 



27r2r 



dk 



{4/3)k^ + dm^/dT 



27r2 



dkk' 



E, 



dT ' 2E, dT 



(24) 



a/, (4/3)fc2 + mf 



dT 



dpo dmi 
dmf dT 



dp 
df 



E 



67r2 



'.4 r 



dfc- 



dT 



fi dm 
2Ef dT 



21 



dpQ dm} 
d^'&T 



In the following, we show results obtained from the 
schematic viscous hydrodynamics with the initial condi- 
tions of To = 0.6 fm/c for the thermalization time, Tq = 
338 MeV for the initial temperature, and Vq = 0.01 c for 
the initial radial flow velocity js^ , which are appropriate 
for heavy-ion collisions at the top RHIC energy. In the 
top panel of Fig. [1] the time evolution of the tempera- 
ture T of the fire-cylinder is given by the solid line. It is 
seen that the QGP, mixed, and HG phases last for about 
5.6 fm/c, 1.5 fm/c, and 10.2 fm/c, respectively, which 
are similar to those in the ideal hydrodynamics with zero 
viscosities in both QGP and HG as shown by the dashed 
line in the top panel of Fig. [TJ For the time evolutions 
of the radial flow velocity Vr — R and transverse radius 
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FIG. 1: Time evolutions of temperature T (top panel), ra- 
dial flow velocity Vr (middle panel), and transverse radius 7? 
(bottom panel) of the fire-cyclinder in viscous (solid lines) and 
ideal (dashed lines) hydrodynamics. 



R of the fire-cylinder in viscous hydrodynamics, they are 
shown by solid lines in the middle and bottom panels of 
Fig. [TJ respectively, and both are slightly above those in 
ideal hydrodynamics (dashed lines). 

Fig. HJa) shows the change of entropy per unit rapidity 
in time. Because of non-zero viscosity, the entropy per 
unit rapidity (solid line) increase with time in the vis- 
cous hydrodynamics, reaching a value at thermal freeze 
out which is significantly larger than that in the ideal 
hydrodynamics (dashed line). Fig. ^h) shows the time 
evolution of the shear tensor components tt^ (solid line) 
and tt'J (dashed line). The two are seen to satisfy ap- 
proximately the relation tt^ « ~T:]tj2. Using this result 
in TT^ -I- TT^ = — tt'''', which is due to traceless of the shear 
tensor, leads to tt^ k, tt''". Since the two are equal in the 
absence of radial flow IM , our results thus indicate that 
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FIG. 2: (a) Entropy per unit rapidity in viscous (solid line) 
and ideal (dashed line) hydrodynamics, (b) Shear tensor com- 
ponents TT^ (solid line) and 7rJ| (dashed line) in the viscous 

hydrodynamics. Dash-dotted and dotted lines are vr^ and Trjj 
in the Navier-Stokes limit. 



For comparisons, we also show in Fig. [IJb) the shear 
tensor components 7r| (dash-dotted line) and Trjj (dotted 
line) obtained in the Navier-Stokes limit, i.e.. 



2vs 



(25) 
(26) 



the assumption of uniform tt^ and tt^ in central heavy- ion 
collisions is reasonable. 



It is seen that they only differ appreciably from those in 
the Israel- Stewart causal viscous hydrodynamics at early 
times. 

Although the evolution dynamics in the viscous hy- 
drodynamics does not differ very much from that in the 
ideal hydrodynamics, the transverse momentum spectra 
of particles can differ significantly. In Fig. [3J we show the 
quark (top panel) and pion (lower panel) transverse mo- 
mentum spectra from the viscous (solid lines) and ideal 
(dashed lines) hydrodynamics at temperatures of 300 and 
150 MeV, respectively. It is seen that the transverse mo- 
mentum spectra from the viscous hydrodynamics are en- 
hanced at high transverse momenta compared to those 
from the ideal hydrodynamics. The inverse slope param- 
eters of the quark and pion transverse momentum spectra 
are 330 and 160 MeV, respectively, which are significantly 
larger than those in the ideal hydrodynamics. The total 
particle density is, however, not much affected by the 
viscosity. To understand these results, we note that the 
particle transverse momentum spectrum in the viscous 
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FIG. 3: Quark (top panel) and pion (lower panel) transverse 
momentum spectra from the viscous (solid lines) and ideal 
(dashed lines) hydrodynamics at T=300 MeV and 150 MeV, 
respectively. 



hydrodynamics is given by [3£ 



/(fc) = /o(fc)+5/(fc) = /o(fc) 



2T{e + p) 



(27) 



where /o(fc) is the equilibrium thermal distribution of 
particles in the ideal hydrodynamics and Sf{k) is the 
viscous correction. In the (r, r, 0, 77) coordinate system, 
the factor PfiPvir^^ can be written as [36| . 

P/jPi/Tt^'^ = m|.cosh^(?; - r;)7r^^ + cos^ {4>p - (fyn''^ 



hp^ sin^ 



•p " <t')T^l + '^T sinh (y - ?7)7r,'J 
2mTPT cosh(?/ — r/) cos((/)p — '?!>)7r'^'', (28) 



where niT = •\/ + with being the transverse 
momentum, and y and (/)p are the energy-momentum ra- 
pidity and the azimuthal angle of the momentum, respec- 
tively. In the absence of radial flow, Eq. pS]) reduces to 



P^iPyT^'^'" = Pt'^'^ 



-pI^I 



(29) 



Since tt'"'' is positive, the particle transverse momentum 
distribution is thus enhanced at high pT- On the other 
hand, the particle number density is not much affected 
by the viscosity as all components of the shear tensor 
have similar magnitude due to its traceless property. 



VI. DILEPTON PRODUCTION IN THE 
VISCOUS HYDRODYNAMICS 

The dilepton production rate from the scattering of 
two particles in hot dense matter is given by 



dN 



(27r)3 (2^)3 /(fci)/(fc^K 



(30) 



where ki and k2 are momenta of the two particles; Urei 
is their relative velocity; and a is the cross section for 
dilepton production from their scattering. For QGP, we 
consider the dominant quark-antiquark annihilation pro- 
cess for dilepton production, and its cross section is 



E 



1 + 2ml /M 



M2 



(31) 



where mq is the quark mass, is the charge of quark 
species i, and ~ (^1 + ^2)^ is the squared invariant 
mass. For HG, the dominant pion-pion annihilation is 
considered, and the cross section is 



A-KO? \F{M 



2\\2 



1 



4m2 



3 M2 V M2 ' 
with the electromagnetic form factor of pion 



■ , ,, v'^i 



M2)2 



2p2 ' 



(32) 



(33) 



where p, p', and p" denote p(770), p'(1450), and 
p"(1700), respectively with their respective width of 



=237 MeV, and F, 



=235 MeV and 
i.802 X 10-^ and 



Vp =153 McV, Vp, 
respective strength N „ = 1, Np 
Npn = 5.93 X 10-3 [if. 

Expressing the relative velocity as Wrci = jki/i?! — 
k2/-E2|, Eq. (|30p can be written in Lorcntz- invariant 
form as 



dN 
d^ 



1 



{2^f 



d^ki d^'ki 
El E2 



/(fcl)/(fc2 



M2 



4m| 
M2 



a(M2), 



(34) 



where uii is the mass of colliding particles. Changing 
variables into = kip^ + k2^ and = (fci^ — kip) 12, we 
can rewrite the Lorentz-invariant phase space as 

d^ki d3k2 d^Vd^'k. d^P Ed^k 



El E2 



E1E2 



E E1E2 



(35) 



with E = El + E2. Because d'^P / E is Lorentz invariant, 
Ed^k/ {E1E2) also should be Lorentz invariant. Boost- 
ing to the center-of-mass frame of the two particles, we 
obtain 



d^P Ed^k rf3p £;'rf3i^/ ^3p 4^31^/ 



E E1E2 E E[E'i E M 



(36) 
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where k' is the momentum, and E[ , £"2 and E' are ener- 
gies of particles 1 and 2, and their total energy in their 
center-of-mass frame, respectively. 

Using Eq.(|27|). the product of particle distribution 
functions becomes 



/(fcl)/(fe) = /o(fcl)/o(fc2)^ 

+ Mk2)Sfih) 



/o(fci)<5/(fc2) 
-Sfih)Sfik2) 



51526 



-P-u/T 



"'2 ""2 " ^^v 



2T2(e+p) 2T2(e+p) 

iLA» 1.1/ L(T l,A 

"-l "-l 'i'2 "-2 ''/ji/''(tA 



4T4(e+p)2 



(37) 



where gi and 32 are degeneracies of particle 1 and 2, 
respectively. Substituting Eqs. ([55)1 and ([57)) into Eq. 
(p4)) . we obtain 



25152 



(2^)6./ ~B 
k''' k"' k''" k 



kf^'k^'lT 



A 



4T4(e + p)2 



A/a 1 



4m| 
IP 



T2(e+p) 
.a(Af2), (38) 
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FIG. 4: (color online) Transverse momentum spectra of dilep- 
tons of various invariant masses in viscous (solid lines) and 
ideal (dashed lines) hydrodynamics. Dotted lines correspond 
to results based on only the first-order correction from the 
modified particle transverse momentum distributions in the 
viscous hydrodynamics. 



where fc^ — {kQ,k'), fc^ = (fcg,— A;') and tt^^, are 
momenta of particles 1 and 2, and the shear tensor 
in their center of mass frame, respectively. By using 
d'^x = TdTrdrdridcf) and d^P/E = irdydP^ and integrat- 
ing over the solid angle of k' in the right hand side of 
Eq. (|38l). the differential yield of dileptons in heavy- ion 
collisions is then given by 



dN 



51.92 



dydM'^dP^ 4(27r 



drr drr 



+ 



M\l + X/S) 
4r2(e-fp) 



/o(a)i^o(/3) 



ij"=0,l,2 



2,3,4 



(39) 



In the above, X = I- Amf/M"^, and Ii{a) and Kj{l3) 
are modified Bessel functions with a = Pt sinh. p/T and 
/? = My cosh p/T, where p = tanh^^ Vr. The coefRcients 
Aij , Bij and details on the derivation of above expression 
are given in the Appendix. 

In Fig. 131 we show by solid lines the transverse mo- 
mentum (pt) spectrum of dileptons from the viscous hy- 
drodynamics for various dilepton invariant masses. Also 
shown by dotted lines are those including only the first- 
order correction from the modified particle transverse 
momentum distributions, i.e., the second term in Eg. pop . 
and they are seen to be very close to those including also 
the second-order correction given by the third term in 



Ea. (|39)) . Compared with those from the ideal hydrody- 
namics given by the first term of Eq. p9p and shown by 
dashed lines, the dilepton spectra in the viscous hydrody- 
namics are enhanced at high pT as in Ref. [IJ] . The vis- 
cous effect on the dilepton transverse momentum spectra 
is thus similar to that on the particle transverse momen- 
tum spectra as a result of enhanced density of quarks 
in QGP or pions in HG at high px- Since dileptons of 
small invariant masses are mainly produced from pion- 
pion annihilations in HG, in which the viscosity is large, 
the viscous effect is thus particularly large for dileptons 
of small invariant masses. Too large a viscous correction 
makes, however, the results from the viscous hydrody- 
namics unreliable. In fact, negative values can appear 
for the dilepton yield from pion-pion annihilation in inte- 
grating Eq. ((39)) over the radial variable r if the dilepton 
invariant mass is large. In this case, we have set the neg- 
ative value to zero. This has, however, a negligible effect 
on the final result as the contribution of pion-pion an- 
nihilations in HG to dileptons of large invariant masses 
is insignificant compared to that from quark-antiquark 
annihilations in QGP. 

Figure O shows the dilepton invariant mass spectra 
from the viscous (solid line) and the ideal (dashed line) 
hydrodynamics. The two are seen to be almost identical 
at all invariant masses as the dilepton transverse momen- 
tum spectra at a fixed invariant mass are similar at low 
transverse momenta, which dominate the dilepton yield, 
in both the viscous and the ideal hydrodynamics. 

In Ref. [1^, it has been shown in the ideal hydrody- 
namics that the differential yield of thermal dileptons 
with a fixed value of transverse mass Mt is essentially 
independent of their transverse momenta if they are pro- 
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of dileptons from QGP at high Pt- 

VII. DISCUSSIONS AND SUMMARY 
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FIG. 5: Dilepton spectra from viscous (solid line) and ideal 
hydrodynamics (dashed line). 



0.01 




1.0 1.5 2.0 
Pj (GeV) 



3.0 



FIG. 6: Transverse momentum spectra of dileptons from QGP 
with transverse mass A'It =1, 1-5, 2, 2.5 and 3 GeV from top 
to bottom. Dashed and solid lines are, respectively, from ideal 
and viscous hydrodynamics. 



duced from QGP [13| , although this so-called AfT-scaling 
is violated for dileptons produced from HG. In Fig. [HI we 
show by dashed lines the transverse momentum spectra 
of dileptons with Mr =1, 1.5, 2, 2.5 and 3 GeV from 
QGP in the ideal hydrodynamics. It is seen that the Mt 
scaling still holds approximately, although we have used 
massive quarks and gluons in our study while massless 
ones were used in Rcf. |13j . The results from the viscous 



hydrodynamics are shown by solid lines and it shows that 
including viscosity leads to a violation of the Mt scaling 



We have derived a set of schematic equations from 
the causal viscous hydrodynamics of Israel-Stewart for 
central relativistic heavy-ion collisions by assuming that 
not only the energy density, pressure and entropy den- 
sity are uniform in the produced fire-cylinder but also 
the azimuthal and space-time rapidity components of the 
shear tensor. Solving these equations using the massive 
quasi-particle model for the equation of state of QGP 
and the resonance gas model for that of HG, we have 
found that the shear viscosity slightly delays the cool- 
ing of produced hot matter and enhances somewhat its 
transverse expansion. It also increases significantly the 
particle distributions at high px, compared with those 
in the ideal hydrodynamics. Using this model, we have 
investigated thermal dilepton production in relativistic 
heavy-ion collisions by including contributions from the 
dominant quark-antiquark and pion-pion annihilations. 
Because of the viscous effect, the dilepton pT spectrum 
is enhanced at high px, which is similar to those found 
in Ref. [l^l based on the non-causal Navier-Stokes vis- 
cous hydrodynamics and the first-order viscous correc- 
tion from the modified particle transverse momentum 
distributions. For the invariant mass spectrum of dilep- 
tons, it is found to differ very little from that in the 
ideal hydrodynamics. We have also studied the effect 
of viscosity on the Mt scaling of thermal dileptons from 
QGP, i.e., the yield is independent of pt for fixed dilepton 
transverse mass Mt, that has been previously predicted 
in the ideal hydrodynamics with massless quarks and glu- 
ons, and it is found that the Mt scaling still holds in the 
ideal hydrodynamics even if QGP is composed of massive 
quarks and gluons as a result of their strong couplings. 
The Mt scaling of dileptons is, however, broken in the 
viscous hydrodynamics due to the enhancement in the 
number density of quarks and antiquarks at high pT by 
the viscous effect. 



Acknowledgements 

This work was supported in part by the U.S. National 
Science Foundation under Grant No. PIIY-0758115 and 
the Welch Foundation under Grant No. A-1358. 



Appendix 

In this Appendix, we give the details on the derivation 
from Eq. (|38|) to Keeping the non- vanishing terms 

in the integration with respect to k' in the right hand 
side of Eq. (l38|) leads to 
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Splitting the integral d'^k' into the radial and angular parts according to 
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and integrating with respect to the solid angle dCl, we obtain 
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where X = I — Am'f/lVP. Using the traceless property of the shear tensor tt^ = 0, Eq. can be rewritten as 
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(43) 



The above integrals can be evaluated by transforming 

c TT^^ to the 
through the 



the shear tensor in the center of mass frame tt^^ to the 
shear tensor in the fire-cylinder frame tt^i^ 
Lorentz transformation, i.e., 
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dy 

ar 



= (a^o coshjy + 0^3 sinh77)(ai,o cosh 77 + 01^3 sinh77)7r^^ 
+ (0^0 sinh?/ + a^3 cosh77)(a^o sinh 77 + a^^ cosh?/)T^7rj,^ 
+(a^i cos (f> + a^2 sin 0)(a^i cos + a^2 sin0)7rrr 
+ (a^i sine/) - a^2 cos0)(a^i sin0 - a^2 cos0)r^7r<^^ 
+ [(a^o cosh ?7 + a^3 sinh 77) (a^i cos cj) + a^2 sin (j)) 
+ {a„o cosh 77 + 0^3 sinh77)(a^i cos0 + a^2 sin (/))]7r7-r, 

(44) 

where a^i/ are components of the Lorentz transformation 
matrix, 
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022 
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(45) 



with 77, y, (/) and being the space-time rapidity, 
energy-momentum rapidity, azimuthal angle in configu- 
ration space and in momentum space, respectively. It is 



straightforward to show that the square brackets in Eq. 
(|43p is a function of y — 77 and 4> — (f>p and leads to mod- 
ified Besscl functions after integration over 77 and 4'p as 
given by Eg.lfM]). where the coefhcients Aij and Bij are 
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with X = Mt/M. 
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